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Abstract 

We present an analysis of the role of the running coupling constant at the intersection of per- 
turbative and nonperturbative QCD. Although the approaches that have been considered so far in 
these two regimes appear to be complementary to each other, a unified description might derive 
through the definition of the effective coupling, as they both provide ways of analyzing its freezing 
at low values of the scale. We extract the effective coupling from all available experimental data 
on the unpolarized structure function of the proton, , at large values of Bjorken x, including the 
resonance region. We suggest that parton-hadron duality observed in this region can be explained 
if nonperturbative effects are included in the coupling constant. The outcome of our analysis is a 
smooth transition from perturbative to nonperturbative QCD physics, embodied in the running of 
the coupling constant at intermediate scales. 
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1. Hard processes are described in QCD by envisaging a perturbative stage (PQCD) where 
a hard collision involving quark and gluons occurs, followed by a nonperturbative stage 
characterizing hadron structure. For example, in Deep Inelastic Scattering (DIS) the hard 
scattering part of the process, 7*g — > X, can be presently described using splitting functions 
up to next-to-next-to- leading order (NNLO) [H, 0], and Wilson coefficients functions up to 
N 3 LO [jjj]. The large distance part is described by the Parton Distribution Functions (PDFs), 
which have to be extracted from fits of high energy experimental data. PDF parametrizations 
along with their uncertainties have been obtained applying this framework up to NNLO, 
by a number of collaborations (see review in jij). NLO PDFs are also important for the 
comparison with data from, e.g., collider processes which are only known up to NLO, and 
when using Monte Carlo event generators. 

Although DIS stands out as the most accurately described process in a QCD factorized 
scenario, at large Bjorken x large logarithms appear in the PQCD description that need 
to be resummed in order to ensure that the perturbative approach can be consistently 
extended. In this region power corrections of both kinematical (Target Mass Corrections, 
TMCs) and dynamical (Higher Twists, HT) origin also affect the extraction of PDFs. They 
are intertwined with Large x Resummation (LxR), and the factorization framework is no 
longer expected to apply straightforwardly jBf. 

Early experimental observations indeed suggested that in specific kinematical regimes 
both the perturbative and nonperturbative stages arise almost ubiquitously, in the sense 
that the nonperturbative description tracks, in average, the behavior of the perturbative 
one. In inclusive eP scattering at large x, a duality between the low-energy and high-energy 
behaviors of a same observable, i.e., the unpolarized structure function, F% was observed 
by Bloom and Gilman who established a connection between the structure function in the 
nucleon resonance region and that in the deep inelastic continuum 0,0|. F%, when averaged 
over the resonance region, was found to be equivalent to the continuation of the deep inelastic 
curve in this region. This concept is known as parton-hadron duality: the resonances are 
not a separate entity but they are an intrinsic part of the scaling behavior of the structure 
function. Note that here duality implies an equality in average between resonances and 
scaling curve, the average being taken at fixed four- momentum transfer, Q 2 > 1 GeV 2 , over 
the same range in the scaling variable, x. This behavior can be taken as signaling a natural 
continuation of the perturbative to the nonperturbative representation. 

Although Bloom-Gilman duality was observed at the inception of QCD, quantitative 
analyses could be attempted only more recently, having at disposal the extensive, high 
precision data from Jefferson Lab (8|. 

In this paper, to evaluate the impact of LxR, we perform an analysis of inclusive eP 
scattering data at large x, including resonance data. To handle the latter, we follow closely 
the analysis of Ref. [9|, where the implications of parton-hadron duality were explored in 
the context of a systematic PQCD based analysis including besides LxR, TMCs and HT 
contributions or, more generally, the evidence for nonperturbative inserts, which are required 
to achieve a fully quantitative fit. Here the relevant kinematical variables are: x = Q 2 /2Mu 
(M being the proton mass and v the energy transfer in the lab system), the four- momentum 
transfer, Q 2 , and the invariant mass for the proton, P, and virtual photon, q, system, 
W 2 = (P + q) 2 (W 2 = Q 2 (l/x - 1) + M 2 ). For large values of Bjorken x > 0.5, and Q 2 
in the multi-GeV 2 region, one has W 2 < 5 GeV 2 , i.e., the cross section is dominated by 
resonance formation. 

As first pointed out in Refs. 0, [Toj, the extension of PQCD evolution to large x must also 
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include LxR effects. A consequence of LxR is that as x increases there is a shift to lower 



values of the scale at which a, is calculated (see for instance Ref. 11 or the classical review 



in Ref. 12j). At large x and low W 2 , this shift requires a freezing of the coupling constant 
in the infrared region. 

The key result in the present work is that, by turning this argument around, or by using 
the fact that our analysis, through LxR, is regulated by the value of the QCD coupling in 
the infrared region, we can extract such coupling from experimental data. In our analysis 
canonical higher-twist terms are suppressed in that their separate contributions to the per- 
turbative curve to be compared to the large x data is negligible. However, nonperturbative 
effects are present as they become absorbed in the coupling's infrared behavior. 

The extracted coupling we obtain is consistent with schemes of scale fixing in which 
a s can be extended to the entire Q 2 domain. Various frameworks have been proposed 



where the effective coupling is free of the Landau pole, e.g., the BLM scheme [13| and its 
recent extension using the Principle of Maximum Comformality [13], methods based on the 
analyticity properties of a s 15, 16], and finally, including nonperturbative effects 17-201 



(see also Refs. [21|, [22[ where a s was defined introducing nonperturbative effects fixed by 
a physical set of parameters). While the focus of the present manuscript is on suggesting 
a way of extracting the running coupling from data, more detailed future studies will be 
dedicated to connecting our approach to the mentioned schemes. 

2. In order to evaluate the effect of LxR we perform a fit of all available large x, eP 
inclusive scattering data. We start from standard parametrizations of the PDFs, and we 
consider systematically the effects of TMCs, and perturbative evolution using either NLO or 
next-to- leading log (NLL) resummed coefficient functions, i.e., with and without LxR. Note 
that the scope of the present fit is not towards a global analysis, but to assess the possible 
interplay among the different components that impact Q 2 evolution at large x, including 
LxR, TMCs, and HTs. In the resonance region, W 2 < 4 GeV 2 , we consider averages of both 
data and theoretical evaluations by comparing limited intervals defined as, 

p m ax(H/ 2 =1.2GeV 2 ) , poxp/ n 2\ 
Dox P /th rn 2N Jz mi n(^ 2 =4GeV 2 ) " XJ 2 ) 

W) rmax (H/ 2 =1.2GeV 2 ) , , Q2) ■ W 

Jzmin(VK 2 =4GeV 2 ) axr 2 l x >V ) 

In the present analysis, we use, for -F 2 cxp , the data from JLab (Hall C, E94110) (23| 
reanalyzed (binning in Q 2 and x) as explained in 24] as well as the SLAC data [IBf. The 



values of Q 2 and the average values of x for each interval are given in Table [H The function 
F 2 th is the theoretical evaluation which is the same in both the DIS and resonance, Eq. ([I]), 
regions. Notice that if Eq. p]) is equal to 1, duality is fulfilled. Since x is integrated over 
the entire resonance region, we are considering global duality. 

The OPE formulation of quark-hadron duality [26j suggests that the higher-twist contri- 
butions to the scaling structure function would either be small or cancel otherwise duality 
would be strongly violated. However, the role of the higher-twist terms is still unclear 
since they would otherwise be expected to dominate the cross section at x — > 1. To an- 
swer the question of the nature of a dual description, two complementary approaches have 
been adopted. The first is the nonperturbative model's view on the scaling of the structure 
functions at low-energies 27-2^]; the second approach consists in a perturbative analy- 



sis 



Id [30|, that through LxR provides a scenario by which the effect of HTs can be 
suppressed in a fully quantitative fit at large x. It is this second approach that we will 
follow in this paper. 
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We evaluate F^ 1 taking into account perturbative evolution at NLO, and introduce sub- 
sequently the effects of TMCs, and LxR. Since only valence quarks distributions are relevant 
in our kinematics, we consider only the Non Singlet (NS) sector, 



F 2 NS ( 




x,Q 2 ) = xq ( x ,Qi) + ^-J2 / dzBUz)- z q{- z ,Q 2 ) , (2) 



The PDFs, q(x,Q 2 ), are taken from current parametrizations. We have chosen to present 
results using the MSTW08 set at NLO as initial parametrization [31?]. We have checked 
that there were no significant discrepancies when using other sets, i.e., CTEQ6 [32] and the 
dynamical GJRFVNS [33]. The function B^ s is the Wilson coefficient function for quark- 
quark. 

By evaluating the ratios i?ex P /th, using current parametrizations, one finds a sensible de- 
viation from the data, even when the theoretical uncertainty from the parametrizations 
is included (Figure [1]). One possible explanation is in the lack of accuracy in the PDF 
parametrizations in the large x, low W 2 domain, since most groups implement much larger 
thresholds for W 2 . The way to a fully quantitative fit would then start from re-fitting the 
large x data with new appropriate sets of PDFs, and simultaneously accounting for both 
TMCs, and LxR. The number of parameters, and the uncertainty associated with this pro- 
cedure would however be dauntingly increasing. For this reason, it is therefore necessary 
to take the preparatory step, conducted with the present analysis, of assessing the relative 
weight of the different contributions. 

The additional corrections due to the finite mass of the initial nucleon, or the TMCs, are 
included directly in F^ s ' th as 34 1 ( se e also review in (HI), 



where F£° = F^ s is the structure function in the absence of TMC. Since TMC should in 
principle be applied also to the HT, we disregard terms of 0(1/Q A ) (36[. Note that the 
expansion in Eq. fl3]) is valid for Q 2 larger than w 1 GeV 2 . TMCs move the ratio closer 
to unity, as rep resented by the open green diamonds in Fig. [TJ Uncertainties on TMCs are 
very small [36]. However a larger error might arise from the procedure used to account for 



TMCs [37|. Studies of the sensitivity to this procedure are on their way and will be published 
elsewhere. At this stage, by including only TMCs and standard PDF parametrizations, we 
still observe a large discrepancy with the data. 

Next, we consider LxR effects. A major consequence of LxR is a shift of the scale at 



which a s is calculated to lower values, with increasing z (see for instance Refs. jllj, |12|, |38| ) . 
This introduces a model dependence within the PQCD approach in that the value of the 
QCD running coupling in the infrared region is regulated by LxR so as to satisfy duality. 
In other words, LxR contains an additional degree of freedom, gathered in the definition of 
the coupling constant, to tune the scaling structure functions. 

LxR arises formally from terms containing powers of ln(l — z), z being the longitudinal 
variable in the evolution equations, that are present in the Wilson coefficient functions 
B^ s (z), in Eq. (J2]). To NLO and in the MS scheme, the Wilson coefficient function for 
quarks reads, 



(4) 
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TABLE I: Upper block: Integrals of JLab data from Refs. 23|, |24j], appearing in the numerator of 
Eq. ([1]). The first column shows the average values of x for each bin. Lower block: SLAC data 251 ] . 
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FIG. 1: Ratio R c * p/th (Q 2 ) of Eq. ([T]) where the theoretical analysis includes PQCD evolution using 
the MSTW08 PDF set (black triangle), and MSTW08 PDF set plus TMCs (open green diamonds). 
Ratios with error bars on data integrated in quadrature ; right panel: ratios with weighted error 
bars on data integrated in quadrature. 



where E.P. means end points and 



Pc 



(0), 



qg 



z) is the LO splitting function for quark-quark. The logarithmic terms, i.e.. 



denotes the standard plus-prescription. The function 

ln(l-z), 

in B^ s (z) become very large at large x values. They need to be resummed to all orders in 
a s . Resummation was first introduced by linking this issue to the definition of the correct 
kinematical variable that determines the phase space for the radiation of gluons at large x. 
This was found to be W 2 



Q 2 (l — z)/z, instead of Q 
of the strong coupling constant becomes ^-dependent 




As a result, the argument 



a s (Q 2 



a. 



'l-z) 



(5) 



In this procedure, however, an ambiguity is introduced, relatedjx) the need of continuing 
the value of a s for low values of its argument, 



i.e.. 



for z — >■ 1 41 



Since the size of this 

ambiguity is of the same order as the higher-twist corrections, it has been considered, in a 
previous work 42], as a source of theoretical error or higher order effects. We investigate 
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the effect induced by changing the argument of a s on the behavior of the ln(l — z)-terms in 
the convolution Eq. fl2]), and resum those terms as 



i r Q2 

ln(l - z) = — - / d\nQ 2 [a s , LO (Q 2 (l - z)) - « s , LO (Q 2 )] = ln LxR 



(6) 



including the complete z dependence of a s>h0 (W 2 ) to all logarithms. 1 Note that we are using 
three different concepts of order expansions. The present analysis is conducted to next-to- 
leading order (we evolve the PDF sets to NLO), to leading-twist (we consider the LT PDFs 
only) and to all logarithms (we include a SjLO (scale) to all logarithms). This resummation is 
easily understood when considering the first term of the expansion of a s (W 2 ) in hx((l—z)/z), 



a s (W 2 ) = a s (Q 2 ) - ^ In \—- J a 2 s (Q 2 ), (7) 
as proposed in Ref. [40|. To all logarithms, the convolution becomes 

F 2 NS ^-(x,Q 2 ) = xq(x,Q 2 ) + £ V f dzBi r -{z) X -q( X -,Q 2 ) , (8) 



where, 



Using ^ V5 ' Hcsum p] us TMCs, in Eq. will make the ratio R decreases substantially, 
essentially leaving no space for HT terms. This is due in our approach mostly to the change 
of the argument of the running coupling constant. At fixed Q 2 , in the integration over 
x < z < 1, the scale W 2 = Q 2 (l — z)/z is shifted and can reach low values, where the 
running of the coupling constant starts blowing up. At this stage, our analysis requires 
nonperturbative information. A way to address this issue is to set a maximum value for 
the longitudinal momentum fraction, z max , which defines a limit from which nonperturbative 
effects have to be accounted for, and to cut a s at the corresponding scale, W 2 (z ma-X ) = 
Q 2 (l — 2 max )/2; max . Larger values of z max correspond to lower values at which the scale should 
be cut in the analysis, meaning that the perturbative value can be used. As we show later, 
large z max occurs in the data at large Q 2 , therefore the effect of the shift in scale gets smaller. 

The functional form ln LxR is therefore slightly changed. Two distinct regions can be 
studied: the "running" behavior in x < z < z m&K and the "steady" behavior 2 max < z < 1, 



F 2 NS ^(x,z m ^Q 2 ) = xq ( x ,Qi) + ^-J2\ / dz [i&CO - Ml"*)' 

+ / dzP^(z) ln LxR + ln LxR , max / dzP^(z) l^g^Q 5 



(10) 



The terms proportional to lnz are not divergent at z — > 1. 
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z mM appears therefore as a free parameter in our analysis. A possible criterion to constrain 
it is to fit the large x data assuming a null direct contribution to the structure function from 
the dynamical HTs namely, for each Q 2 bin we define z max by varying _R ox p/ th as a function 
of z max , so that 



i? cxp/th (^ ax ,Q 2 



dxFf p (x,Q 2 ) 



dx F< 



NS.Ke 



Jcxp 



(x, 2 max , Q 



(11) 



In Eq. (FIT]) . F 2 W5 ' Rcsum (a;, z max , Q 2 ) was evaluated including TMCs, resummation to all log, 
and setting possible dynamical HT contributions to zero. The latter get, however, absorbed 
in the coupling's infrared behavior. More precisely, the suppression of HTs in the structure 
function is compensated by the behavior of a s in the infrared region. As a result, contrarily 
to what originally deduced in, e.g., Ref. j43|, a definite role of nonperturbative corrections 
is obtained, pointing at the fact that duality, defined on the basis of a dominance of single 
parton scattering, i.e., suppression of final state interactions, might indeed be broken. 

Results are represented by the red hexagons in Fig. [2J The integrals values are given in 
Tab. HT1 together with the corresponding values for z max . Since, for the largest values of Q 2 , 
Q 2 = 10, 15GeV 2 , outside the resonance region, on Fig. [fl z max becomes closer to 1, we do 
not consider those data points in what follows. 
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FIG. 2: The ratio i? cxp/th (x ave , Q 2 ) as a function of Q 2 , on the left pannel, and as a function of x avc 
on the right. Same as Fig. [1] but with the red hexagon representing the LxR results of Tab. HU 
The key shows the results corresponding to JLab data. The open triangle, full diamonds and open 
hexagons corresponds to SLAC data. 



3. Based on the results of our analysis of large x data including TMCs and LxR, we now 
extract a s by assuming that it runs from the onset of a minimal scale which is determined 
from the comparison with data, and it is frozen from that minimal scale downward to the 
real photon limit (scale=0 GeV 2 ). As one can see from Table HT] data in the resonance region 
are crucial for this determination. 

In Fig. E]we show our extracted value a SjNLO (scale) where we used the MS scheme outside 
the IR region, for the same value of A throughout this paper. a s was obtained as an exact 



solution to NLO [2l|. Our theoretical error band is defined by the shift in z max from the 
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TABLE II: Integrals at each stage. 
/ Rcsum (z max ,Q 2 ). 



In the last columns: the value z r , 



associated with 



different bins displayed in Table [TT] namely. 

k l ^max,t) 
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S,NLO 



Qt 
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1, ... 10 



(12) 



i corresponds to the data points. Including this error band, our extracted frozen value of 
the coupling constant is, using the MSTW08 PDF set for the analysis, 



0.1337 < 



ft, 



, (scale ->■ OGeV 2 



< 0.1839 



7T 



(13) 



In the figure we also report values from the extraction using polarized eP scattering data 



in Ref. [44j-|47|. These values represent the first extraction of an effective coupling in the 
IR region that was obtained by analyzing the data relevant for the study of the GDH sum 
rule. To extract the coupling constant, the MS expression of the Bjorken sum rule up to the 
5th order in alpha (calculated in the MS scheme) was used. In order to compare with our 
extraction using the F$ observable, the finite value for a s (0) found in (45h47] was rescaled in 
44 1 assuming the validity of the commensurate scale relations (47[ in the entire range of the 
scale entering the analysis. The agreement with our analysis which is totally independent, 
is impressive. 

4. In conclusion, we presented an extraction of a s using eP scattering data at large x. A 
careful analysis of all the contributions appearing at large x including TMCs and LxR, was 
performed. The central value for a s {Q 2 < lGeV 2 )/-7r was found to be 0.1588. This value is 



46 



in agreement with the extraction from the GDH sum rule analysis \AA 

When considering PQCD observables at low scales, we implicitly face an interpretation 
problem. In the multi-GeV 2 region and at large-x, where the resonances lie, perturbative 
QCD is pushed to its limits. Both higher terms in the perturbative expansion of that 
observable, and power corrections need be taken into account. In the present approach, this 
transition is taken into account by re-interpreting the running coupling constant, at the scale 
of transition instead. By tuning the scaling structure functions to the averaged data in the 
resonance region, we parametrize the realization of duality through an infrared fixed-point 
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scale (GeV 2 ) 



FIG. 3: (color online) Extraction of a s . The blue dashed curve represents the exact NLO solution 
for the running coupling constant in MS scheme. The solid blue curve represents the coupling 
constant obtained from our analysis using inclusive electron scattering data at large x. Owing 
to large x resummation, at lower values of the scale, a s = a SjN Lo (min) is frozen as explained 
in the text. The grey area represents the region where the freezing occurs for JLab data, while 
the hatched area corresponds the freezing region determined from SLAC data. This error band 
represents the theoretical uncertainty in our analysis. We also plot results extrapolated 44] from 
the recent analysis in Refs. 45l. l4q] : the red squares correspond to a s extracted from Hall B CLAS 
EGlb, with statistical uncertainties; the orange triangles corresponds to Hall A E94010 / CLAS 
EGla data, the uncertainty here contains both statistics and systematics. 



below which the strong coupling constant stops running. The value we have found for the 
transition scale is situated around 1 GeV 2 . The uncertainties on that value are, however, 
still important. 

We report an interesting observation that by connecting the values of the coupling using 
two different observables and schemes, namely the g\ scheme from the GDH sum rule extrac- 
tion, and our MS based extraction, we obtain an excellent agreement by simply extending 



49]). While our conclusion 



the commensurate scale relations to the IR region (see also [48|, 
ensues from a perturbative analysis, in the near future, we will explore the role of nonper- 
turbative effective couplings as well. An interesting avenue is provided in this direction by 
the evaluation of the strong coupling constant in the infrared region using the AdS/CFT 
conjecture that was recently performed in Ref. 47 . 



The importance of finite couplings has been highlighted many times in the Literature, 
see, e.g., [14J and references therein. Our analysis, would allow in principle, to extract from 
a fit the nonperturbative parameters often present in the proposed functional forms for the 



running of a s , e.g., in Refs. |15|, [17|, [18 
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